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SHIELDING APPLICATION OF PERTURBATION THEORY TO DETERMINE CHANGES 
IN NEUTRON AND GAMMA DOSES DUE TO CHANGES IN SHIELD LAYERS 

by Daniel Fieno 
Lewis Research Center 

SUMMARY 

Perturbation theory formulas were derived and applied to determine changes in 
neutron and gamma-ray doses due to changes in various radiation shield layers for fixed 
sources. For a given source and detector position the perturbation method enables dose 
derivatives due to all layer changes to be determined from one forward and one inhomog- 
eneous adjoint calculation. Hence, the perturbation method for obtaining dose deriva- 
tives requires fewer computations for design studies of multilayer shields. 

An illustrative problem is presented in which a comparison was made of the frac- 
tional change in the detector dose per unit change in the thickness of each shield layer as 
calculated by perturbation theory and by successive direct calculation. The problem in- 
volved a neutron source surrounded by spherical layers of tungsten and lithium hydride. 
The perturbation theory gave -0. 238 cm”'*' for the fractional change in the neutron dose 
per unit change in the tungsten layer thickness; the result from direct calculations was 
-0. 234 cm"*. Similar close agreement was obtained for the effect of a change in the 
thickness of the lithium hydride layer on the neutron dose and for the effects of changes 
in both the tungsten and lithium hydride layers on the detector dose due to secondary 
gamma rays generated in the tungsten layer. 

INTRODUCTION 

Nuclear reactors for space propulsion systems or auxiliary power systems require 
shielding to protect personnel, system components, and/or payload from neutrons and 
gamma-ray radiation. These shields are usually composed of alternate layers of heavy 
metal and hydrogeneous material. One shielding problem for these systems is to deter- 
mine the number, arrangement, and thickness of these shield layers that will result in 
minimum-weight shields while maintaining the radiation field at various positions around 
the system within prescribed limits. 
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One -dimensional shield weight optimization procedures have been evolved for the 
design of layered configurations (refs. 1 to 3) and these have been extended to constraints 
in several directions around two-dimensional configurations (refs. 3 to 5). All these re- 
quire the determination of the change in dose at a detector due to a change in thickness of 
each layer of a given shield configuration. These are generally obtained with S n trans- 
port calculations, requiring at least as many separate calculations as there are shield 
layers. As the optimization proceeds and the configuration changes, this procedure has 
to be repeated until the configuration converges to an optimum one. The optimization 
procedure becomes time consuming and, particularly for two-dimensional configurations, 
would require excessive amounts of computer time. 

The present work applies perturbation theory for systems with fixed sources to ob- 
tain information to apply to multilayer -shielding problems. The use of perturbation 
theory in reactor physics applications is not new (ref. 6), the primary usage being the 
determination of reactivity effects. The application of perturbation methods to nonreac- 
tivity problems appears mainly in the Russian literature (refs. 7 and 8). The perturba- 
tion method utilizes only one forward and one inhomogeneous adjoint solution for an 
initially selected layered shield configuration to determine all necessary dose deriva- 
tives at a detector for a given neutron or gamma-ray source component. The method 
appreciably reduces the amount of calculations required for multilayer shields. In addi- 
tion, perturbation theory yields a considerable amount of other information concerning 
the effects of individual reaction processes (scattering, absorption, etc. ) on the dose 
rate. 

The equations for this shield perturbation method are developed in the appendix and 
are applied to a simple spherical system to determine the dose derivatives of the various 
layers with respect to the layer thickness. This illustrative problem also serves to 
demonstrate the ability of the perturbation method to extract additional informative data 
concerning interaction processes of the various shield materials. 


METHOD OF ANALYSIS 

This section presents the equations used in the perturbation analysis for systems 
with fixed sources. The theory is outlined for a fixed- source neutron system since the 
results for fixed-source gamma systems and coupled neutron-gamma systems are simi- 
lar in form. A more complete derivation of the equations is given in the appendix. The 
symbols used are mainly conventional; those that are not are defined. 

Consider the time-independent multigroup form of the inhomogeneous neutron trans- 
port equation. This equation is given by 
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fi- V0 g (r, fi) + CT g (r)0 g (r, fi) = Qg(r, S2) 


Z 2 Z _+_1 V' 

477 / J 

1=0 g- 


l 

CT Z (r,g'-g) • ^ 

m=-Z 


<m^ P Z,m^ 


g = 1, 2, 


( 1 ) 


where d> (r, f2) is the angular flux for group g, o is the macroscopic total cross sec- 
tion for group g, Q^(r, ft) is the angular fixed source for group g, cr^g'-g) is the l - 
order expansion coefficient of the scattering from group g' to group g, (pf (r) is the 
(l, m) 11 expansion coefficient of the angular flux for group g in a spherical harmonics 

E i- U 

series, and P^ m (ft) is the (Z, m) n spherical harmonic as a function of the unit direc- 
tion vector ft. The boundary condition for the angular flux at the convex boundary of 
the system is taken as 


<t> g (r s ,fi") = ° g = 1, 2, . . G (2) 

This boundary condition is the usual condition and implies that no neutrons enter the 
system from the vacuum. Here r„ indicates the outer boundary of the system and 0“ 
designates an inward direction. 

Once the solution, 0 g (r, ft), to equation (1) has been obtained, then a physical 
process which depends linearly on the angular flux may be written as a functional, F. 

The general multigroup form of F is given by 

0 g (r, ft)P g (r, ft) drdft (3) 

g fi H 

The functional F gives contributions caused by an interaction process for the neutrons 
in the system. The quantity P g (r, G) has the role of an interaction cross section for 
group g at position r for neutrons moving in the direction ft. The integrations are 
over all directions and over the entire system designated by R. 

If a functional at some specific point, r d , is required, equation (3) may be written 
as 
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(4b) 


F(r d ) = 



S2)P g ( r d> fi) 


where 6(r - r^) is the Dirac delta function. The general expression, equation (3), will 
be used in subsequent derivations. Equation (4a) will be invoked when results such as 
radiation dose at a given position are required. 

Following Usachev (ref. 7) an equation can be derived for the importance of neutrons 
for the physical process described by the functional F. The importance i// (r, fi), which 

o 

is equal to the effect F at some detector position or region due to a unit source of 
neutrons placed in group g at (r, f2), is given in multigroup form by 


-S3- Vi// g (r, fi) + CT g (r)^ g (r, to) - P g (r, to) 


i=0 g' m =-l 


m 


(fi) g=l, 2, . 


(5) 


The boundary condition for the importance xf/Ar, fi) follows from the fact that, if a neu- 

& 

tron escapes from the system, it cannot contribute to the physical process of interest. 
Therefore, its importance must be zero and the boundary condition becomes 


* g (r s , ^ + ) = 0 g = 1, 2, . . . , G (6) 

where fi + designates an outward direction. 

The next step in the analysis is to derive an equation in terms of a perturbation in 
the density of the material constituents of the system. The perturbation in the density 
of a material, which can be related to a change in the thickness of that material shield- 
ing layer, is introduced in the equations through the macroscopic cross sections and ap- 
pears as a perturbation in the cross sections of the system. If equation (1) were solved 
directly with the perturbed cross sections, a perturbed value O) for the angular 

flux would be obtained. The transport equation for the perturbed system is 


ft) + jvgfr) + 6CT g (r)J U) = Q^(r, H) + 6Q^(r, i?) 

l 

g’-g) + 6c/(r, g’-gjj • g = 1, 2, . . . , G (7) 

l =0 g' m=-£ 
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along with the boundary condition 


^(r g , 0“) = 0 g = 1, 2, . . . , G (8) 

In the preceding equation, .^(r, fi) represents the perturbed value of the angular flux for 
group g. Perturbations in the cross sections are indicated by 6a. To generalize the 
analysis, the source term Qg(r, fi) is assumed to be perturbed by 6Q g (r, £2). 

The neutron transport equation for the perturbed system (eq. (7)), along with the 
equation for the importance (eq. (5)) and the boundary condition (eqs. (6) and (8)), can be 
manipulated algebraically to give the following result: 




Z 2 1 + 1 

4tt 

1 = 0 




dr i// g ( r, fi)Q g (r, O) - 


y J1 dfi f 


dr i// g (r, fi)6Q^(r, ft) 


+ 



dr J^(r, fi)P g (r, ft) = 0 


(9) 


Equation (9) is an exact relationship and represents a condition which the perturbed 
angular flux and importance must satisfy. If all the perturbations in equation (9) are 
zero, then equation (9) along with equation (3) gives the important result 



The perturbed flux .^(r, S2) has been replaced by its unperturbed value 0 g (r, fi). The 
result expressed by equations (3) and (10) is well known (ref. 9) and is a necessary con- 
dition which the angular fluxes and importances must satisfy. 

In the perturbed system the functional F p for the physical process of interest is 
equal to 


F p = y C f rlf, fi)P p (r, a) drdf2 (11) 

“ •'n y R s g 

g 

Here the general case of a perturbed system is considered for which P p (r, fi) is the per- 
turbed value of the quantity P (r, Q). The quantity P p (r, fi) is taken to be 

o o 

P P (r, O) = P g (r, O) + 6Pg(r, a) (12) 

Next consider the difference 6F of the functionals F p and F. Using equations (10), 
(11), and (12) in equation (9) results in 

6F = F p - F 



Equation (13) is the desired perturbation expression for a system with fixed sources. 
This equation gives the change in the functional F for a given physical process due to 
changes in various neutron cross sections, changes in the fixed source, and changes in 
the physical process itself. This equation can also be used to give the change in F due 
to perturbations in the individual material constituents of a shield layer or due to per- 
turbations in individual neutron cross sections. Equation (13) is an exact relationship 
since no approximations are inherent in its derivation. However, since only low -order 
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perturbation theory is being considered, the perturbed flux .^(r, O) may be replaced 
by its unperturbed value $g(r, £2). 

In shielding analysis, generally the biological dose rate at a detector position r^ 
is the functional of greatest interest. The interaction cross sections Pg(r, £2) are the 
flux-to-dose conversion factors which are independent of r and £2. Thus, we can ex- 
press P (r, £2) as 

o 


d 

P Jr, £2) = S- (14) 

s 4n 

where d^ is the flux-to-dose conversion factor for group g. Insertion of equation (14) 
into equation (4b) gives 


F ’ j; £ <> 

g 

Thus, the dose functional F is 1/4 tt times the dose rate at position r^. 

The first two terms on the right side of equation (13), which depend on the density, 
can be put in terms of a material efficiency function h q (r) (ref. 8). That is, these 
terms can be defined as 


47r6F CT 


£ f h q (r)6p q (r) 
q ^ 


dr" 


(16) 


where q indicates a given material component and SF^ is the contribution of these 
terms to QF. The efficiency function h q (r) defines the change in the functional F 
caused by the introduction of a unit quantity of material q at position r. The efficiency 
functions are derived in the appendix for the various terms which depend on the neutron 
cross sections. The efficiency functions are useful for determining the effect on F of 
the different material components of a region at any position within that region; also 
they are useful for determining the effect on F of the important neutron interactions 
which occur for a given material. 

The calculational procedure used in applying the perturbation method to determine 
changes in neutron and gamma doses due to changes in shield layers is as follows: 

(1) For an initial, layered shield configuration, solve equation (1) along with the 

boundary condition (eq. (2)) for the angular flux </> (r, £5). 

6 


7 



(2) Solve equation (5) along with the boundary condition (eq. (6)) for the angular dose 
importance function i£/g(r, ^). This solution is for an inhomogeneous adjoint source at a 
detector location and for a set of flux-to-dose conversion factors. 

(3) Compute the dose functional F by equations (3) and (10). The functionals F, ob- 
tained from these two equations, must be equal; differences are caused by the numerical 
approximations used in obtaining the angular fluxes and dose importances. 

(4) Compute, by using equation (13), the change in the dose functional 6F caused by 
perturbations in the cross sections (due to density perturbations), perturbations in the 
fixed source, and perturbations in the dose response functions. The quantity 6F is 
computed separately for each shield layer for a nominal change in the layer density 5p. 

(5) The effect of each shield layer on F can be expressed as the fractional change 
in the dose functional per unit change in layer thickness 6F/F6t. The coefficient 
5F/F5t is determined by taking 6F, computed for a nominal change in the density of a 
shield layer 6 p, and dividing by the product of the functional F and the change in the 
thickness of the shield layer 6t. This change in the thickness of a shield layer 6t is 
related to a change in the density of a shield layer by 6t = Spt/p where p is the density 
of the shield layer and t is its thickness; that is, by taking the product pt as a con- 
stant. 


ILLUSTRATIVE PROBLEM 

Some of the capabilities of the perturbation method are illustrated on the one- 
dimensional spherical system shown in figure 1. This system has a central cavity or 
void region with a radius of 30 centimeters followed by a 15-centimeter-thick layer of 

tungsten (W) and a 60-centimeter-thick layer of lithium hydride (LiH). A neutron source 
17 

of 1. 7x10 neutrons per second having a fission spectrum was uniformly distributed in 
the cavity. Doses and dose importance functions were evaluated for a detector located 
0. 25 centimeter from the outer boundary of the system. The composition of the system 
is listed in table I. 

Neutron fluxes and importance values throughout the system were determined by 
using the S n method to solve both the forward and appropriate inhomogeneous adjoint 
equations. A 26 -energy-group split was used to represent the neutron spectra. Elastic 
scattering was treated through the Pg order. An S^g Gauss- Legendre quadrature was 
utilized in these calculations. 

The perturbation method was applied to secondary gamma sources - treating sources 
from neutron captures and from inelastic scatters individually. This necessitated gen- 
erating each of these sources from neutron absorption and neutron inelastic scatter re- 
action rates in the tungsten layer (secondary gamma sources in the lithium hydride layer 
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Figure 1. - One-dimensional spherical system used in illustrative problem. 
(Detector position at 0. 25 cm from outer LiH boundary. ) 


TABLE I. - COMPOSITION OF 


SPHERICAL SYSTEM 


Region 

Element 

Density, 

3 

g cm 

Atom density, 
atom b-cm 

1 

(a) 



2 

W 

19. 27 

0. 0631 


w 182 

5. 04 

0. 0167 


W'83 

2. 77 

. 0091 


W 184 

5. 92 

.0194 


w 186 

5. 54 

.0179 

3 

LiH 

0. 748 

0. 0568 


Li 6 

0. 042 

0. 0042 


Li 7 

. 612 

. 0526 


H 

. 094 

. 0568 


a Void. 



were neglected) and evaluating separately capture gamma and inelastic scatter gamma 
fluxes; the importance function for the gamma dose is independent of the gamma source 
component and is evaluated only once. S n transport calculations using a 15-energy- 
group split to represent the gamma spectra were made in a manner similar to those 
made for neutrons. In general, the perturbation method could more easily be applied 
by using a 41 -energy -group coupled neutron-gamma S n transport calculation if the in- 
tent were to determine effects of shield layers on the combined neutron and gamma dose 
rates. 

The effect of varying the density of each layer can be expressed as the fractional 
change in dose functional per unit change in layer thickness 6F/F6t. These quantities 
are evaluated and compared with results obtained from successive direct calculations 
in which nominal changes in the thickness of each layer were made. The effects of 
various interactions in both shield layers on these dose derivatives are evaluated, as 
well as the influence of various material components of each layer. 

Additional information is obtained from evaluating the perturbation equations for the 
neutron and gamma dose functions. This information is given in terms of a dose effi- 
ciency function which is the derivative of the dose functional with respect to weight and 
is evaluated as a function of position in each layer. 


RESULTS 

Effect of Tungsten and Lithium Hydride Layers on Neutron Dose 

9 

The neutron dose rate was calculated from equation (3) to be 2. 90x10 rem per 

17 

hour for the source of 1. 7x10 neutrons per second. The neutron dose calculated from 

o 

equation (10) was computed to be 3. 11x10 rem per hour. The ratio of the dose calcu- 
lated from equations (10) and (3) was equal to 1. 07. This ratio should be unity, and 
departure from unity is a measure of the approximations inherent in the S n method for 
a given spatial mesh and angular quadrature. Thus, if the values of the dose functionals 
differ by a large enough amount to preclude obtaining the accuracy desired, then S n 
calculated angular fluxes and importances must be recalculated with a finer spatial 
mesh and angular quadrature. 

For the illustrative problem considered, the importance function (r) gives the 
dose at the detector due to a unit isotropic spherical shell source of neutrons located at 
position r in energy group g. Since the qualitative features of this neutron dose im- 
portance function \F (r) are less well-known that those for the flux $ (r), the neutron 

o o 

importance function obtained as a solution to equation (5) is shown in figure 2 for groups 
1, 5, 10, and 15. Results for the other neutron energy groups show similar spatial 
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Dose importance function for neutrons, 5glr), rem/hr/neutron in group g 



Radius, cm 

Figure 2. - Dose importance function for neutrons $ q (r), giving 
response of unit isotropic spherical shell sources. IDose de- 
tector at 0. 25 cm from outer boundary. ) 


11 






characteristics. Figure 2 indicates the strong effect of the LiH on the variation of the 
neutron dose importance function and the marked effect with energy group. 

The effect of the tungsten layer on the neutron dose is presented in table II in terms 
of a coefficient 5F/FSt which gives the fractional change in the neutron dose per cen- 
timeter. The data in table II give the effect of the important reactions which a neutron 
can undergo with nuclei of tungsten. Inelastic neutron scattering accounts for 61 percent 
of the contribution to the effect of the tungsten layer, elastic neutron scattering contrib- 
utes 20 percent, (n, 2n) reactions contribute 19 percent, and neutron absorptions have a 
negligible contribution. 


TABLE II. - NEUTRON DOSE DERIVATIVES 
6F/F6t FOR TUNGSTEN LAYER, 
COMPUTED BY PERTURBATION 


THEORY 


Neutron reaction 

Fractional change in 
neutron dose per 
centimeter, 
6F/F6t, 
cm"* 

Elastic scattering 

-0. 04 7 

Inelastic scattering 

-. 146 

(n, 2n) Reaction 

-.045 

Absorption 

-. 0001 

Total 

-0. 238 


The effect of the lithium hydride layer on the neutron dose is presented in table III 
in terms of the coefficient 6F/F5t. In addition, partial effects due to constituents hy- 
drogen (H) and lithium (Li) are also presented. The data in table III indicate the dom- 
inant role of elastic neutron scattering for low -mass-number elements on the neutron 
dose functional. The data in the table also show that the effects of H and Li are nearly 
the same with respect to the neutron dose. 

A comparison is given in table IV of the fractional change in the neutron dose func- 
tional per unit change in the tungsten layer thickness SF/F6t between the perturbation 
theory result and the result calculated by successive direct transport calculations. The 
perturbation theory result is -0. 238 cm"\ while the result from direct flux calculations 
is -0. 234 cm'l The two results are in close agreement with each other. For the LiH 
layer the perturbation theory result for the coefficient 6F/FSt is -0. 138 cm"\ while 
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TABLE III. - NEUTRON DOSE DERIVATIVES 


6F/F6t FOR LITHIUM HYDRIDE LAYER, 


COMPUTED BY PERTURBATION THEORY 


Neutron reaction 

Material 

H 

Li 

LiH 

Fractional change in neu- 
tron dose per centi- 
meter, 6F/F6t, cm’* 

Elastic scattering 

-0.070 

-0.040 

-0. 110 

Inelastic scattering 


-. 027 

-.027 

Absorption 


-.0005 

-.0005 

Total 

-0. 070 

-0. 068 

-0. 138 


TABLE IV. - COMPARISON - FOR THE NEUTRON 
PROBLEM - OF DERIVATIVE 6F/F6t 
CALCULATED BY PERTURBATION 
THEORY AND BY SUCCESSIVE 
DIRECT TRANSPORT 
CALCULATIONS 


Region 

Fractional change in neutron 

Ratio of pertur- 


dose per unit change in 

bation result to 


layer thickness, 6F/F6t, 
cm’ 1 

direct result 


Perturbation 

Direct calcu- 



calculation 

lation 



■>wm- 


1. 02 




. 99 


the result obtained from direct flux calculations is -0. 140 cm'l Again the two results 
are in close agreement with each other. 

In computing the perturbation theory results, material efficiency functions h(r) 
were generated as intermediate results for each material component and for each neu- 
tron reaction as a function of position within a layer. This additional information can be 
used in other aspects of shield analysis since h(r) is a derivative giving the change in 
the dose functional F per unit mass of material introduced at position r. For this il- 
lustrative problem the total h(r) of each material component is shown in figure 3(a) for 
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£» 7 

the W of the W layer and in figures 3(b) and (e) for the H, Li , and Li components of 
the LiH layer. 


Effect of Tungsten and Lithium Hydride Layers on Secondary Gamma Dose 

The gamma-ray dose rate was calculated from equation (3) to be 2. 59x10^ rads per 

hour for a gamma-ray source due to neutron absorptions in W of 3. 98x10 photons per 

second. By using equation (10) the gamma dose was also calculated to be 2. 59x10^ rads 

per hour. The gamma-ray dose was calculated from equation (3) to be 7. 10x10** rads 

1 7 

per hour for a gamma-ray source due to neutron inelastic scattering in W of 8. 17x10 
photons per second. From equation (10) the photon dose for this case was calculated to 
be 7. 19x10** rads per hour. 

The gamma-ray dose importance function ^gC?) is shown in figure 4 for energy 
groups 1, 5, and 10. Results for the other photon energy groups show similar spatial 
characteristics. Figure 4 indicates the strong effect of the W on the variation of the 
gamma-ray importance function. 

The secondary source distributions due to W(n, y) and W(n, n'y) reactions are shown 
in figure 5(a). The perturbation formula (eq. (13)) for the secondary gamma source in- 
cludes a term which is based on the difference between the perturbed and unperturbed 
neutron fluxes. This perturbed secondary source term is estimated in the following 
manner: The neutron transport problem is calculated as a two-outer-iteration problem. 
The first iteration is performed for a small perturbation in the density of the W layer. 
This yields a set of perturbed neutron fluxes which are stored. The density of W is next 
restored to its unperturbed value, and the calculation continues as a second outer itera- 
tion for the unperturbed neutron fluxes. The difference between the perturbed and un- 
perturbed values of the neutron fluxes are then used to calculate the perturbations in the 
secondary sources. Figure 5(b) presents the perturbations in the secondary source dis- 
tributions for the W(n, y) and W(n, n'y) reactions for a 1-percent increase in the density 
of the W layer. In a coupled neutron-gamma problem the secondary gamma sources 
never explicitly appear in the equations and the preceding source perturbation calcula- 
tion is not necessary. 

The effect of the tungsten layer on the gamma doses is presented in table V in terms 
of the coefficient 5F/F6t. The data in table V give the effect of each reaction which a 
gamma ray can undergo with a nuclei of W. In addition, the effect on the gamma dose 
of perturbations in the secondary gamma-ray source is also given. The data show that 
these perturbations in the secondary gamma source are a major contribution to the 
gamma dose functional. The effect of the LiH layer on the gamma doses is presented 
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Secondary photon source, gammas/cnr-sec 



(a) Secondary photon sources in tungsten layer (b) Perturbation in secondary gamma sources 

caused by neutron (n,y) and (n, n'y) reactions caused by a 1-percent increase in tungsten 

with tungsten nuclei. density. 

Figure 5. - Secondary source distributions. 


TABLE V. - SECONDARY GAMMA DOSE DERIVATIVES 
6F/Ffit FOR TUNGSTEN LAYER, COMPUTED BY 
PERTURBATION THEORY 



Secondary photon source 

W(n,y) reaction 

W(n, n'y) reaction 

Fractional change in photon dose 
per centimeter, 6F/F5t, cm'* 

Photon reaction: 



Scattering 

-0. 041 

-0. 068 

Absorption 

-. 029 

-. 022 

Secondary source 

-.040 

-. 082 

perturbations 



Total 

-0. 110 

-0. 172 




TABLE VI. - SECONDARY GAMMA DOSE DERIVATIVES 6F/F6t FOR 


LITHIUM HYDRIDE LAYER, COMPUTED BY 
PERTURBATION THEORY 



Secondary photon source 


W(n,y) reaction 

W(n 

n'y) reaction 


Material 


H 

Li 

LiH 

H 

Li 

LiH 


Fractional change in photon dose per unit change in 
layer thickness, 6F/F6t, cm"' 

Photon reaction: 
Scattering 
Absorption 

-0.0065 

(a) 

-0. 0195 
(a) 

-0.026 

(a) 

-0. 010 
(a) 

-0.029 

(a) 

-0. 039 
(a) 

Total 

-0.0065 

-0.0195 

-0. 026 

-0.010 

-0.029 

-0.039 


a Negligible. 


TABLE Vn. - COMPARISON - FOR GAMMA 
PROBLEM - OF DERIVATIVE 6F/F6t 
CALCULATED BY PERTURBATION 
THEORY AND BY SUCCESSIVE 
DIRECT TRANSPORT 
CALCULATIONS 


Region 

Fractional change in photon 
dose per unit change in layer 
thickness, 5F/F6t, cm"' 

Ratio of pertur- 
bation result to 
direct result 

Perturbation 

calculation 

Direct calcu- 
lation 

Secondary photon source: W(n, y) reaction 

W 

■ 



LiH 


1 

msM 

Secondary photon source: W(n, n 

’y) reaction 

W 

-0. 172 

-0. 168 

1. 02 

LiH 

-.039 

-. 039 

1. 00 








in table VI in terms of the coefficient 5F/F6t for the two secondary gamma-ray 
sources. 

The coefficient 6F/F6t as computed by perturbation theory and successive direct 
flux calculations is compared in table VII for the W and LiH layers. The table shows 
that the coefficient 6F/FSt computed by perturbation theory is in agreement with the 
result obtained from successive direct flux calculations for all cases. 

The total material efficiency functions h(r) with respect to the gamma-ray dose are 
shown in figures 6(a) and (b). These functions were obtained by using the perturbation 
equations for each material component of the LiH layer and for the W of the W layer. 



-lOxlO" 2 



(a) Tungsten. Valuesof hfrl are positive from (b) Lithium hydride. 

30- to 34-centimeter radius and of order 10 . 

Figure 6. - Total photon dose efficiency function for tungsten and for lithium hydrode. Secondary gamma source: W(n.y) reactions. 
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CONCLUSIONS 


The result for the coefficient 6F/F6t giving the fractional change in the dose per 
unit change in shield layer thickness that were obtained by using perturbation theory 
formulas for fixed sources were in excellent agreement with those obtained by using 
successive direct transport calculations. The perturbation method also gives a con- 
siderable amount of information concerning each reaction which a neutron or photon 
can undergo with the nuclei of various materials and thus adds to the qualitative under- 
standing of changes in shield design. This same information is not easily or readily ob- 
tainable from successive direct transport calculations. Although the features of this 
perturbation method were exhibited with a simple one-dimensional spherical problem, 
the extension of the perturbation method to two-dimensional geometry and to coupled 
neutron - gamma-ray problems poses no additional difficulty. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 1, 1972, 

503-05. 
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APPENDIX - PERTURBATION THEORY FOR SYSTEMS WITH FIXED SOURCES 

Neutron Transport Equation 

The application of perturbation theory in the computation of reactivity effects is well 
known and, in fact, was used from the very beginning of the atomic energy program. 

Less well known is the application of perturbation theory to the computation of nonreac- 
tivity effects such as, for example, the breeding ratio of a reactor system, the sensi- 
tivity of the biological dose rate to variations in shield composition, and the sensitivity 
of other reactor processes to variations in the system. The theory underlying the ap- 
plication of perturbation methods to such nonreactivity problems appear mainly in the 
Russian literature (refs. 7 and 8). Lewins reviews, in reference 6, this and other de- 
velopments in perturbation theory. This appendix gives the details of the mathematical 
analysis underlying the use of perturbation theory for systems with fixed sources. The 
symbols used are mainly conventional; those that are not are defined. 

To provide a theoretical basis for understanding perturbation theory for systems 
with fixed sources, consider the time -independent multigroup form of the neutron trans- 
port equation. This equation is given by 


ft- V0 g (r, ft) + CT g(r)0 g (r, ft) = Q g (r, ft) + — ^ <x( r, g'-g)<f> g ,(r) 

4 “ g’ 

+ X ? 17'X s ‘ iS<? ’ g, " g) ' X ^' m(?)P '- m(S> 8=1,2 G <A1> 

l-l g' m=-Z 

where <£ (r, ft) is the angular flux for group g, cr is the macroscopic total cross sec- 

O — ► — O £ T'h 

tion for group g, Qg(r, ft) is the angular fixed source for group g, o (g'~g) is the l 
order macroscopic cross section for scattering from group g' to group g, cjjf (r) is 

f'h ^ ? *** 

the {l , m) 1 expansion coefficient of the angular flux of group g in a spherical harmonic 
. series, and (ft) is the {l , m) n spherical harmonic as a function of the unit direc- 

tion vector ft. 

For convenience in performing subsequent manipulations, the equation is written for 
one material component. In reality, all the terms involving cross sections in equa- 
tion (Al) are summations for all the material components present at position r. Al- 
though equation (Al) will be studied from the point of view of neutron transport, it is ap- 
parent that the analysis will also pertain to photon transport or coupled neutron and 
gamma-ray transport when a proper identification of terms is made. This is essential, 
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for example, in determining the biological dose rate external to a reactor shield which 
is composed of both a neutron and photon component. 

In deriving equation (Al) the angular flux which occurs in the elastic scattering 
integral on the right side of the equation has been expanded in a series of spherical 
harmonics and is given by 


0 g (r, fi) = 



2 rhS ft 


An 


*F f m< r > p «,m< 0) 


(A2) 


l =0 m=-Z 


The spherical harmonics are defined here as 




t 


( 2 ' 6 mo )(Z " m) 


,-, 1/2 


(l + m)! 




• P ™(p) cos(m <p) l > 0, 0 < m < Z (A3a) 


and 


P i m (0) - P (Z ~ m ^ - l 1/2 • P™(p)sin(m <p) 

L ’ m l{l + m)! J L 


l > 0, 0 < m < £ 


(A3b) 


where P“ (p) are the associated Legendre polynomials of argument p and where 6 mo 
is the Kronecker delta which has the property that 


J kk' 


k = k’ 


k * k* 


(A4) 


For the local Cartesian coordinate system at position r shown in the following 
sketch (where 6 is the polar angle and cp is the azimuthal angle), p = cos 0. 


Z 
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Thus, from the sketch, the unit vector fl is given by 

£2 = (A5) 

and the element of solid angle d£2 given by 

dfi = d/idcp (A6) 

For geometries with symmetry in the azimuthal angle cp, all the terms in sin(m(p) are 
identically zero in the expansion of the angular flux. Thus, the complexity of equa- 
tion (A2) is reduced for simple geometries. For example, only the terms in Q (f2) - 
the Legendre polynomials - are needed for one -dimensional slab or spherical geometry. 

The spherical harmonics defined here have two important properties. One of these 
is the orthogonality property given by 




(A7) 


The other property is the addition theorem, which is given by the equation 


p,<5-5-)=£ 

m=-Z 


(A8) 


By the orthogonality property of these spherical harmonics, the m (r) of equa- 
tion (A2) can be determined. These are given by 


■V ? ’ aP *, m (5)dn 


(A9) 


The quantity q(t) is the scalar neutron flux for group g and will be represented more 
simply as (r). The net neutron current vector J (r) is given by the equation 

s & 


J (r) = L 4> (r, n)£2 d£2 
& * / n & 


= x (r) + j _j(r) + k Q (r) 


S 


(A10) 
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where i, j, and k are unit vectors along the x-, y-, and z-coordinate axes, respec- 
tively. . 

The unit vector is given in terms of the spherical harmonics by 

S2=lP 1 + _ 1 (fi) + kP lj0 (fi) (All) 

Equations (A10) and (All) can next be combined to give the dot product of 0 and J g (r). 
This is 


J g^ = Z 0 1 m ( " )P 1, m (fi) (A12) 

m=-l 

To complete the description of equation (Al), we must recall that a number of quan- 
tities must satisfy certain normalization requirements. These are 


es, g 


-z 

g' 


S“ s ( g- g ') 


(A13a) 


S in, g = Z (A13W 

g' 

S n2n, g = | Z S n2n< g - g ’) <A13c> 

g’ 

( ^l) g 2 es, g = Z 2 es^~^ ^ (A13d) 

nr * 


< V l < g L» g ' S es, g = Z S e S < g - g ’> <A13e) 

g’ 

Note also that, when neutron leakage is treated by means of a buckling term, the quantity 

ct (r) must be suitably modified. The quantities E . Z- . and Z ~ are the ma- 
5 g iii, g nzn, g 

croscopic elastic, inelastic, and (n, 2n) cross sections for group g, respectively. The 

quantities { /i L ) and <P^(/j. L )) represent the average cosine and the average l th 
g g 

order Legendre polynomial of elastic scattering for group g, respectively. 
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Functional for a Physical Process 


A physical process which can occur in a system described by equation (Al) and which 
depends linearly on the angular flux can be written as a functional F which is given in 
multigroup form by 


<f> g ( r, s5)P g (r, Si) drdft (A14) 

g g R 

where F is given by 

o 

F = / / <p (r, ft)Pjr, ft) drdft (A15) 

& */R ® & 

The functional F gives contributions caused by the physical process for the neutrons in 
the system. The quantity P (r, S2) has the role of an interaction cross section for the 
physical process occurring in group g at position r for neutrons moving in the direc- 
tion ft. The integrations in equations (A14) and (A15) are over all directions ft and 
for the entire system designated as R. 

If a functional at some specific point, r d , is required, equation (A14) can be written 
as 


F(F d ) - LAX 0 g ( r, ft)P g (r, ft)6(r - r d ) drdft (A16a) 

g O R 

F(? d ) = <P g (r d , fi)P g (r d , ft) dft (A16b) 

g ^ n 

where 6(r - r d ) is the Dirac delta function. The general expression, equation (A14), 
will be used in subsequent derivations. Equation (A16a) will be invoked when results 
such as radiation dose at a given position are required. 

In shielding analysis, as an example, generally the biological dose rate at a detec- 
tor position r d is the functional of greatest interest. The interaction cross sections 
P (r, ft) are the flux-to-dose conversion factors which are independent of r and ft. 

Thus, P (r, ft) can be expressed as 

o 
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(A17) 


P„(r, O) = i 

& 471 

where is the flux-to-dose conversion factor for group g. Insertion of equation (A17) 
into equation (A16b) gives 


F = 


_L 

4 77 


Z 

g 


(A18) 


Thus, for this shielding application, the functional F is 1/477 times the dose rate at 
position r^. 


Importance Equation for a Given Physical Process 

Following Usachev (ref. 7), an equation can be derived for the importance of neu- 
trons for a given physical process. The importance of an initial distribution of neutrons 
is equal to the sum of the importance of all the scattered neutrons and of the number of 
processes which occur. Since the importance of each process is, by definition, unity, 
the preceding statement results in a law of conservation of importance. Therefore, let 
i//g(r, Q) be the importance of a neutron introduced in group g at position r and moving 
in the direction fi. If r' corresponds to the position of a neutron in group g after it 
travels a distance ds in the direction at position r, then the following expression 
gives the importance due to M neutrons: 


Mip (r, 57) = M\p (r',57) ■ (Probability that a neutron does not have a collision in ds) 

& & 

+ (Importance due to scattered neutrons produced by original neutrons in distance ds) 


+ (Importance ol physical process in distance ds) 


(A 19) 


The terms on the right side of equation (A19) which are expressed in words can be easily 
given in terms of multigroup cross sections, that is, 
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I 


M^ g (r, S2) = M • 


[i - a g (r)ds]i// g (r', S2) + 



CT (r, g~g')* g ,(r) 


Mds ' Y ' Y ^ g " g,) ' Y ^ ' m ® Pi . “ (S> 

g' m=-Z 


1=1 


+ Mds • P (r, n) • (Unit importance) 

D 


(A20) 


Next expand \p ( r', O) in a Taylor's series about r: 
© 


_ _ 3^ (r, 9 

i// (r', 0) = xp (r, S2) + — a ds + 0(dsr 

g g 3s 


(A21) 


Also make the following replacement: 


3ij/ g (r, O) 
3s 


= n- V^/ g (r, S2) 


(A22) 


Substitute equation (A21) into equation (A20), divide through by Mds, and let ds — 0. 
The result will be an integrodifferential equation for the importance of neutrons to the 
given physical process. This equation is given by 


-fi-Vi// g (r, O) + o- g (r)^g(r, n) = -L ^ a(r, g-*g')^ g ,(r) 


g’ 


OO 


7 =1 


2^+1 
4 77 



l 

m =-l 


+ P g (r, fi) g = 1, 2, . . . , G 


(A23) 


Although equation (A23) is, in the strictest sense, not adjoint to equation (Al), it is 

adjoint-like in that it gives the importance of neutrons for the physical process. Note 

that different interaction processes described by the function P g (r, S) result in different 

importance functions \p (r, S2). 

g 
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The boundary condition for the importance i^g(r, £2) follows from the fact that, if a 
neutron escapes from a system, it cannot contribute to the physical process. Therefore, 
its importance must be zero and the boundary condition becomes 


* g (r s , f2 + ) = 0 g = 1, 2, . . . , G (A24) 

where r indicates the outer boundary of the system and S2 + designates an outward 
direction. 

In deriving equation (A23) and in analogy with equation (Al), the importance function 
which occurs in the elastic scattering integral on the right side of the equation has been 
expanded in a series of spherical harmonics and is given by 


°° l 
2=0 m =-l 


21 + 1 
477 




(A25) 


By the orthogonality of the spherical harmonics used here, the pf (r) of equation (A25) 

C j ill 

can be determined. They are given by 



™ Cl) ‘A ^g (:r ' aP !,m® dn 


(A26) 


The quantity q(t) is the scalar importance for group g and will be represented more 
simply as (r). A vector J+(r), which is the net importance current vector, is de- 

o o 

fined by the equation 


Jg( r ) = ^ g (r, d S2 = fyg j(r) + _](?) + 0 (r) 

u 6 

Equations (All) and (A27) can be combined to yield the dot product of S2 and 
that is, 


(A27) 


(r), 


n- Jg(r) 




m=-l 


(A28) 
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Perturbation Equation for a Given Physical Process 


Consider a system which is perturbed in a manner such that all the cross sections 
are varied. The solution of equation (Al) with the perturbed cross sections yields a 
perturbed value for the angular flux. The neutron transport equation for this case is 
given by 


n- fi) + [og(r) + oa(r)J- J^(r, O) = — • ^ [o<r, g'-g) + 6cr(r, g’-g)J • J|.,(r) 

47 1 g’ 

°° l 

+ g ’~ g) + 6S ^ (? ’ g ’^ g) ] ‘ X ^ 

2=1 g’ m=-l 

+ Qg(r, fi) + SQ g (r, S2) g = 1, 2, . . . , G (A29) 

In the preceding equation $F{ r , fi) represents the perturbed value of the angular flux for 

g 

group g. Perturbations in the cross sections are indicated by replacing a cross section 

a by a + 6cr where 6cr is the perturbation. The source term Q^(r, S2) is assumed to 

undergo a perturbation in the amount of 6Q g (r, Si). 

Multiply equation (A29) by i// (r, SI) and then integrate the resulting equation over 
, & 

all drdfi as well as summing the results over all groups g. Next multiply equa- 
tion (A23) by S ? " (r, Si) and then integrate the resulting equation over all drdfi as well 

o 

as summing the results over all groups g. Finally, subtract the second equation from 
the first; and the result, after considerable algebraic manipulation, will be 




tre 6 V~ r) V r - S) 


_(r , n dS + 


6u(r, s ,- li)'l'g(r) 


°° L 

Z~ ir • Z / dF Z< m<T ' ) ' Z w< ^ mfr) 


L - 1 m=-£ 


z dn X dr *g (r ’ S) ^ (r ’ S) ■ z dn x dr v r> c)5Q K (r ’ n > 

E g 


Y A dn f dr ^(r,5)P (r,jj) =0 
A' •'n At g g 


(A 30) 


29 



In the first term of the preceding equation, dS is an element on the outer surface S of 
the system. 

If all the perturbations in equation (A30) are zero, then the equation reduces to 



0 g (r s , S 2 )if/ (r s , J 2 )S 2 - n ds - 



+ y C dn f dF 0 (r, B)Pjr, fi) = 0 (A31) 

Sr g e 

g 

where JT. (r, fi) has been replaced by its unperturbed value $ (r, fi). Equation (A33) is 
g g 

well known and is the multigroup form of a similar equation as obtained by Hansen and 

Sandmeier (ref. 9). 

The boundary condition for the angular flux at the convex boundary of the system 
will be taken as 


0 g ( r s ,sT) -0 g = 1, 2, . . ., G (A32) 

This boundary condition is the usual condition implying that no neutrons enter the system 
from the vacuum. Here designates an inward direction. Using boundary conditions 
(A24) and (A32) reduces equation (A31) to 


F - 






dr i// g (r, n)Q^(r, fi) 


(A33) 


Equation (A33) shows that the functional F for a given physical process may be 
computed in two different ways. The equation may also be used to give another useful 
description for the importance function i// g (r, fi). Let the source Qg(r, Q,) be given by 



The source given by this equation describes a source of neutrons in group g' at posi- 
tion r Q and moving the direction £2 Q . Next substitute equation (A34) for Qg(r, S2) in 
equation (A3 3). This will give 


d£2 jT dr 0 g (r, n)P g (r, £2) (= F) (A35) 

Equation (A3 5) leads to the conclusion that the importance function i^ g (r, £2) is the func- 
tional for a physical process for a unit source of neutrons placed in group g at posi- 
tion r and traveling in direction £2. In other words, the importance function j^ g (r, £2) 
is a Green's function. Further, a solution of equation (A23) provides the functional for 
a physical process for every possible unit source. Equation (A3 5) gives the units of the 
importance function ^(r) as that of the functional F per source particle. 

From equation (A33) we can deduce the following additional conclusions concerning 
the functional F: If the value of the functional is needed for a number of different source 
distributions Q^(r, 12), then it is more efficient to solve for the importance function 
^/^(r, £2) and perform the convolution indicated by the second part of equation (A33). 
However, if the value of the functional is needed for a number of different interaction 
cross sections P g (r, fi), then it is more efficient to solve for the fluxes <p g (r, £2) and 
perform the convolution indicated by the first part of equation (A33). 

In the perturbed reactor the functional giving the contributions caused by the given 
physical process for the neutrons in the system is, in analogy to equation (A14), equal to 

F p =T f f jr (r, £2)P p (r, £2) drd£2 (A36) 

“ Jo. SR g g 

Here the general case of a perturbed system is considered for which P p (r, £2) is the per- 
turbed value of the cross-section-like quantity P (r, £2). The quantity P p (r, £2) is taken 

& o 

to be 

pP(r, £2) = P g (r, £2) + 6P g (?, S) (A37) 

where 5P g (r, £2) represents the perturbation in the interaction cross sections Pg(r, £2). 

Next consider the difference SF of the functionals F p and F. Eliminate the first 
term of equation (A30) by using the boundary conditions given by equations (A24) and 
(A32). Then insert equations (A33) and (A36) into equation (A30) giving the following 
result: 
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4tt 6F = 4n-F p - 4;rF 


dfi j~ dr 6CTg(r)^/g(r, Q) jr(r, Q) 


00 

Z( 2z + i) 'E L dr 6cr es^ r ’ g '‘* g ^f : m^ r ^ 

l =0 m =-l g’ g 

dr Yj i g’^ ) ' Yj 6CT in^’ g '~ g) ^g® + Yj S g'^ ' X 6<T n2n^’ g'-g^gf?) 


4/7 * dn dr 4/g(r, Q)5Qg(v, Q) + An ■ %A dft J" dr j/r(r, fi)6P g (r, fi) (A38) 


Equation (A38) is the desired perturbation expression for a system with fixed sources. 
This equation gives the change in the functional F for a given physical process as a 
result of changes in various neutron cross sections, changes in the fixed source, and 
changes in the interaction process itself. The third term on the left side of equa- 
tion (A30) has been split into its component parts comprising elastic scattering, inelas- 
tic scattering, and the (n, 2n) reaction. The elastic scattering reaction is considered to 

7 j-i, 

be anisotropic, with Scr" (g'—g) being the perturbed values of the l -order macro- 
scopic transfer cross sections. Inelastic neutron scattering and the (n, 2n) reaction are 
considered to be isotropic in the laboratory coordinate system with Scr^Cg’— g) and 
5CTn2n(§ ,_ *g) being the perturbed values of the macroscopic inelastic scattering and (n, 2n) 
cross sections, respectively. 

Equation (A38) is an exact relationship since no approximations are inherent in its 
derivation. However, since only low-order perturbation theory is being considered, 
the perturbed flux & (r, Q) can be replaced by its unperturbed value, <fi (r, f2). 

6> o 

Efficiency Function for Functional F 

Definition of efficiency function . - Abagyan et al. introduce in reference 8 the con- 
cept of a material efficiency function h q (r) relative to a functional F: 
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h q^x _ 4n-5F q (r) 
6p q (r) 


(A39) 



The efficiency function defines the change in a functional F caused by the introduction 
of a unit quantity of material q at the point r. The efficiency function shows the sen- 
sitivity of the physical process to changes in the density of the different materials in the 
system. As can be seen, the quantity h q (r) has the dimension of the change in the func- 
tional per gram at position r for material q. 

If the efficiency functions h q (r) are known, then the variation of the functional 6F q 
with the change of density of material q in a unit volume at position r by a value 
6p q (r) is equal to 


47r5F q (r) = h q (r)6p q (r) 


(A40) 


and 


4t7SF^ = f' h q (r) 6p q (r) dr 


(A41) 


4tt5F = Yj 47T6F q (A42) 

q 

Thus, a knowledge of efficiency functions provides information about changes in func- 
tionals resulting from small changes in material densities 5p q (r). The efficiency func- 
tions h q (r) can be calculated by means of the perturbation expression given by equa- 
tion (A38). 

Equation (A38) is the basic equation for the computation of the changes in the func- 
tional F due to changes in the cross sections. Some of the terms in equation (A38) can 
be put in terms of efficiency functions as defined by equation (A39), as will be shown in 
the next sections which follow, with the perturbed fluxes replaced by their unperturbed 
values. 

Effi ciency function for total-cross-section term . - Let the change in the functional 
F due to changes in the total cross section for material q (the first term on the right 
side of eq. (A38)) be designated as 6F^. That is, let 

47T6F? = -477- V dS2 f dr 6cr q (r)^/ (r, S2)0 (r, C2) (A43) 

1 v Jsi Jr g g s 

g 

The integration over in equation (A43) can be easily performed by using equations 
(A2) and (A25) along with the equation for the orthogonality property of the spherical 
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harmonics, equation (A7). This will give for 6F^ the following expression: 


4ff6f ? = - £ £ d? 6a g^ • Yj (2i + 1) * Yj *1, m® (A44) 

g K 1=0 m =-l 


The quantity 6 ct^ is the change in the macroscopic total cross section for material 

o O 

q due to changes in the density of material q or to changes in the cross section itself. 
In other words, 5 cr^ is equal to 

b 


* a -~ ’ av g - “2® g 
A q 


(A45) 


where p? is the microscopic total cross section for group g for material compo- 
nent q. The quantity is the atomic weight for material q, N ay is the Avogadro 
number, and N^(r) is the original atom density of material q at position r. The 
quantity Sp^ may be, for example, the estimated experimental error in the total 
cross section for group g for material q. The quantity 6p^(r) represents the change in 
the density of material q at position r. 

Let us write equation (A43) in the following form: 

4ff6Fj = 4ff5F^(p) + 4jr6Fj?(XS) (A46) 

That is, we wish to consider 5F^(p) as the change in the functional SF^ due to a change 
in the density of material q and to consider 6F^(XS) as the change in the functional due 
to changes in the group -dependent microscopic-cross-section data. Thus, 47rSF^(p) is 
given by 



with the efficiency function h^(r) being equal to 

^t,g- Z ( 21 + !)■ X * 1 , m (7) 

1=0 m =-l 


h?(r) = 


N. 


av 


A^ 


E 


(A47) 


(A48) 
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1 


The efficiency function h^(r) is defined here as the change in the functional F caused 
by a unit change in the density of material q at position r for the effect due to the 
total cross section. The component parts of the efficiency function given by equa- 
tion (A48) as a result of absorption, elastic scattering, inelastic scattering, and (n, 2n) 
reactions may be obtained by using the appropriate partial cross section in place of the 
total cross section. 

The quantity SF^(XS) is the change in the functional due to changes in the total 
cross section itself and is given by 


4tt6F^(XS) = 


f drj^Njr) ■ £ 64? g * Z & + « * Z «f, m^Z , m G>] 

K L g Z=0 m =-Z J 


(A49) 


As indicated previously, the effect on the functional F of the partial cross sections 
which make up the total cross section can be computed by using equation (A49) with the 
appropriate cross section in place of the total cross section. 

Efficiency function for elastic scattering transfer terms . - Let the change in the 
functional F due to changes in the elastic scattering transfer cross sections for ma- 
terial q (the second term on the right side of eq. (A38)) be designated by 6F^ . Then 
this quantity is equal to 


°° Z 

W es = Z (2Z + • Z L d? Z *?,' ' Z 6 4s, cf’ m® 

Z=0 g’ g (A50) 


7 4-L 

The quantity 5S1_ Jr, g'— g) is the change in the macroscopic Z -order elastic scat- 
es ? q 

tering transfer cross section for group g’ to group g for material q due to changes in 
the density of material q or to changes in the transfer cross section itself. In other 
words, 


<s, q (? ’ g, ~ g) = ~ ’ 4s, q (g ’~ g) + N o^> 6 4s, q (g, ~ g) (A51) 

A^ 

where u (g 1 — g) is the microscopic Z -order elastic scattering transfer cross sec- 

es, q 7 

tion for group g’ to group g for material component q. The quantity 6ju„ c (g 1 — g) may 

gs ? q 

be associated with, for example, experimental errors in the elastic scattering transfer 
cross sections. 
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Let us write equation (A50) in the following way, 


47r5F q = 4?r6F q (p) + 4?r6F q (XS) 


(A52) 


1 


That is, we wish to consider SF q g (p) as the change in the functional 
change in the density of material q and to consider 6F q (XS) as the 

tVi 6S 

tional due to changes in the l -order microscopic elastic scattering 
sections. Thus, 47r6F q (p) is given by 

C fa 


6F q due to a 

cS 

change in the func- 
transfer cross 


4tt 6F q g (p) = f h q g (r)6p q (r) dr 

•'R 


(A53) 


with the efficiency function h q (r) being defined by 


N 

h q „(r) =-JZ. 


es' 


A q Z=0 


E ® + ' Z Z *f, m® • E 4s, m ® (A54) 


m=-Z g’ 


The efficiency function h q (r) is defined here as the change in the functional F caused 
by a unit change in the density of material q at position r for the effect due to the 
elastic scattering transfer cross sections. 

The quantity <5F q (XS) is the change in the functional due to changes in the elastic 

Gfa 

scattering transfer cross sections and is equal to 




c t> 


g' 

l , m 



g 


6 '4s,q (g ’- g) 

J 

(A55) 


Equations (A50) to (A55) also apply to inelastic scattering and the (n, 2n) reaction 
given by the third and fourth terms on the right side of equation (A38). All that needs 
to be done is to replace the elastic scattering transfer cross sections by appropriate 
inelastic scattering or (n, 2n) reaction transfer cross sections. Although the various 
scattering processes are split for convenience, we could combine all the scattering 
terms (elastic, inelastic, and (n, 2n)), which are second, third, and fourth terms of 
equation (A38). The result is equations (A54) and (A55) with /ii„(g'— g) replaced by 
!± (g'— g) - the total microscopic transfer cross section. 

Efficiency function for the fixed source term . - Let the change in the functional F 
due to changes in the fixed source term (fifth term on the right side of eq. (A38)) be 
designated by 6Fq. Then this quantity is equal to 
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4tt6Fq = 47r^ dS2 dr ^ g (r, S2)6Q g (r, O) 


(A56) 


The quantity 6Q g (r, fi) represents the change in the fixed source for group g. Equa- 
tion (A56) can be evaluated directly if values of 6Qg(r, fi) are known. 

In computing the biological dose rate, if both the neutron and photon components of 
the dose rate must be computed, the following analysis is necessary (this does not apply 
to coupled neutron-gamma problems): Part of the photon dose rate is the result of 
gamma rays from radiative neutron capture and inelastic neutron scattering. This 
source term is assumed to be isotropic and is equal to 


where Q q (r) is given by 


Qg( r) = fgQ^C?) 


(A57) 


G 


'N 


Q q (r) = ^ a q gr(r)4> gI (r) 
g’=l 


(A58) 


In equation (A58), Q q (r) represents the process rate per cubic centimeter for material 
q at position r" for a reaction with cross section a for neutron energy group g. 

X J 5 Q 

The summation extends over all neutron energy groups. The quantity rj is the number 

o 

of photons which appear in gamma group g for each reaction which occurs for mate- 
rial q. 

The quantity SQ^C?) now represents a change in the secondary gamma source Q q (r) 
due to changes in the density of material q and the reaction cross sections n 
Thus, 6Q q (r) is equal to 

o 


X, g 


5Q q ff) = 6(f|Q q (T)} - ■ Q q (r) + f q - Q q (J) 

p q (r) 


G 


+ f 


g 


N 

Z 

g'=l 


a x,g* (7)6 \^) 


+ f 


J N 


q N^C?) 


I 

g’=l 


6m x, g' 




(A59) 


Equation (A59) is obtained by using equation (A58) and an equation analogous to equa- 


In equation (A59), 5f q is the change in the number of photons 

o 


tion (A45) for 6cr q (r). 

■'S & 

which appear in gamma group g for each reaction which occurs for material q. 


The 


37 



quantity 5p q is the change in the microscopic reaction cross section for neutron 
S , 

energy group g for material q. The quantity 6<f>g(r) is the perturbation in the flux and 
is equal to ^(r) - $ (r). 

Let us write equation (A56) for a given material q in the following way: 

4tt6F^ = 4tt6Fq(p) + 4tt6F^(XS) (A60) 

That is, we wish to consider 6F^(p) as the change in the functional 6F^ due to a change 
in the density of material q and to consider 6F^(XS) as the change in the functional due 
to changes in the reaction cross sections and gamma spectra for the particular reaction 
which produces the secondary photon source. Thus, 47r6F^(p) is given by 


4jt6 Fq(p) = jT h^(r)6p q (r) dr 


= f dr<jv f^ (r)l ■ Q<1(?) + V cr ,(r) • t Cr)l 

7r lg J l p(r) g^l J 

with the efficiency function h^(r) being equal to 


(A61) 


= jQ^Cr) . 1 


X N 


h q (7) 

Q 1 P( r) 5p 


g'=i J g 


(A62) 


The efficiency function h^(r) is defined here as the change in the functional F caused 
by a unit change in the density of material q at position r resulting from a change in 
the secondary source Q q (r). 

The quantity <5F^(XS) is the change in the functional due to cross-section changes 
and is given by 


47r6Fi(XS) = / dr 

y •'R 


I 

Lg 


6f q ^ g (r) 


• Q q (r) 


+ 



G 


N 


^ 5 ^x,g'' *g'® 

g' = l 


(A63) 
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Total efficiency function for a material . - The total efficiency function h^(r) for a 
material q is the sum of the efficiency functions derived in the preceding sections for 
various terms of equation (A38), that is, 

h^(r) = h q (r) + h q g (r) + hj? n (r) + h q 2n (r) (A64) 

The efficiency function h^,(r) is, therefore, the total change in the functional F per 
gram of material q introduced at position r. The total dose functional 6F^(p) caused 
by a change in the density of material q is equal to 

47r6F^,(p) = / h q (r)6p q (r) dr (A65) 

•s R 


or equivalently, 


6F q (p) = 8F q (p) + 6F q g (p) + 6F q n (p) + 6F q 2n (p) (A66) 

Perturbation Equation for Functional F in Terms of Various Partial Functionals 

The total change in the functional F for a system with a fixed source is given by 
equation (A38). The preceding analysis has described the use of efficiency functions to 
calculate changes in the functional due to the density changes in the material components 
of the system. The change in the functional F due to changes in the cross sections 
6Fj(XS) of material q is given by 

6F q (XS) = 6F q (XS) + 6F q g (XS) + 6F? n (XS) + 5F q 2n (XS) (A67) 

Finally, equation (A38) can be evaluated in terms of equations (A56), (A66), and 
(A67). The result is 


4tt6F = Y, 4ff 6F q (p) + Y 4?r6F q (XS) + 4 ttSF q + 47r6F p (A68) 

q q 

where 6F p is defined by the equation 
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(A69) 


4jr5Fp = ■ 



The summations in equation (A68) are over all material components q. 


Calculational Procedure 

The calculational procedure for using perturbation theory for systems with fixed 
sources is as follows: 

(1) Solve equation (Al), by using the boundary condition given by equation (A32), for 
the neutron fluxes using the given fixed source terms. 

(2) Solve equation (A23), by using the boundary condition given by equation (A24), 
for the importances for the interaction process of interest. 

(3) Compute the functional F for the interaction process of interest by using equa- 
tion (A14). 

(4) Compute the functional F by using equation (A33). The two values of this func- 
tional, computed by using equations (A14) and (A33), must be equal and represent a 
necessary condition on the solutions obtained for the fluxes and importances. 

(5) Compute the change in the functional 477SF for the physical process of interest 
due to perturbations in the cross sections, perturbations in the fixed source, and per- 
turbations in the interaction process itself by using equation (A38). 
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